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Homodyne detection for atmosphere channels 
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We give a systematic theoretical description of homodyne detection in the case where both the 
signal and the local oscillator pass through the turbulent atmosphere. Imperfect knowledge of the 
local-oscillator amplitude is effectively included in a noisy density operator, leading to postpro- 
cessing noise. Alternatively, we propose a technique with monitored transmission coefficient of the 
atmosphere, which is free of postprocessing noise. 
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I. INTRODUCTION 

Long-distance quantum communication necessar- 
ily deals with strong unwanted effects of the environment. 
In this context, one usually compares two types of chan- 
nels: optical fibers and free space. For purposes of quan- 
tum optics, it is important that in fibers [3] one usually 
deals with a stable intensity attenuation and with strong 
depolarization effects. In free-space channels [3,01, the 
situation is different: The attenuation randomly fluctu- 
ates and the depolarization effect is negligibly small. 

An important method for measuring the quantum-light 
characteristics is the technique of balanced homodyne de- 
tection @, [zl ■ In this case, the signal field is combined 
through a 50:50 beam splitter with a strong coherent 
field, the local oscillator. The difference of photocounts 
in two outputs of the beam splitter is proportional to the 
field quadrature. By applying this procedure for different 
values of the local-oscillator phase, one could get com- 
plete information about the quantum state of the signal. 
Particularly, one can reconstruct the density operator in 
different representations 0-Il|- 

The application of homodyne detection for long- 
distance quantum communications in free-space channels 
meets the problem of phase synchronization between the 
signal and the local oscillator. A possible way to over- 
come this difficulty could be based on the technique of 
the optical frequency comb In this case, the de- 

tected signal will be randomized by the atmosphere with 
respect to both the amplitude and the phase A 
more traditional way to provide such a synchronization 
is to derive the signal and the local oscillator from the 
same source. However, in this case the local oscillator 
will also be affected by the atmospheric turbulence. At 
least part of this problem can be resolved by sending the 
signal and the local oscillator from the same source in or- 
thogonally polarized modes jl^ . Since the atmospheric 
depolarization effects are negligible, the phase synchro- 
nization is not destroyed in such an experiment. On the 



E-mail address: |sem@iop.kiev.ua] 



other hand, in this case the local-oscillator amplitude 
randomly fluctuates due to the atmospheric noise. As 
a result, the problem is how to connect the photocount 
difference with the field quadrature. 

In the present paper we consider the situation when the 
signal and the local oscillator pass through the turbulent 
atmosphere in orthogonally polarized modes. First, we 
analyze the scheme proposed in Ref. In this case, 

the photocount difference can be connected with the field 
quadrature by using a certain reference value of the local- 
oscillator amplitude. This is equivalent to the use of a 
reference transmission coefficient of the atmosphere, for 
example its mean value. This results in a kind of noise, 
which is related to the postprocessing procedure. Re- 
constructed with such a procedure, the density operator 
may even fail to satisfy the fundamental requirement of 
positive semidefiniteness, which is a serious disadvantage 
of this method. To resolve this deficiency, we propose 
a procedure with a permanently monitored transmission 
coefficient. This renders it possible to recover the true 
values of the field quadratures from the measured photo- 
count differences. In this case, the postprocessing noise 
disappears. 

The paper is organized as follows. In Sec.|lTl we derive 
an expression for the statistics of photocount differences 
for the scheme considered in Ref. [12'| . This result is used 
in Sec. IIIIl where the photocount difference is connected 
with the field quadrature by using a fixed reference trans- 
mission coefficient. A method based on monitoring the 
transmission coefficient is developed in Sec. lIVI In SecfVl 
we derive input-output relations for the normally ordered 
covariance matrix and consider the effect of quadrature 
squeezing of the light passing through the atmosphere. 
A summary and conclusions are given in Sec. I VII 



II. STATISTICS OF PHOTOCOUNT 
DIFFERENCES 

Let us consider an experimental scenario as imple- 
mented in Ref. [l2| with the local oscillator copropagat- 
ing with the signal field in the same spatial but differ- 
ent polarization modes; see Fig. [TJ The half-wave plate 
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HWP and the polarization beam-splitter PBS at the re- 
ceiver prepare 50:50 field superposition of the signal and 
the local oscillator as done in the standard homodyne de- 
tection. The detectors Di and D2 are used for measuring 
the photocount difference An, which is used for further 
analysis. 

The probability distribution of the photocount dif- 
ference An in the Heisenberg picture is given by (cf. 
Ref. [13) 



= Tr 



noisy 



(1) 



where p is the input-signal density operator and K^^^ 
is the noisy positive operator- valued measure (POVM) of 
photocount differences. The latter, 
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is expressed in terms of the POVM of photocounts. 



where 



■exp {-rih^^i - Nn 



(3) 



• : denotes the normal ordering prescription. 
The coherent-state vector \re^'^') represents the local os- 
cillator of amplitude r and phase (/?. The vacuum-state 
vector |0) includes all modes of the environment. More- 
over, r] is the detection efficiency, and Nnc is the mean 
value of noise counts caused by stray light as well as dark 
counts; see Ref. [31 ■ The annihilation and creation op- 
erators bi and SJ, respectively, represent the light modes 
at the ith output of the polarization beam splitter PBS. 
This beam splitter and the half wave-plate HWP can be 
described by the input-output relations. 



(4) 
(5) 



where 5s and bio are annihilation operators of the signal 
and the local oscillator at the input of the beam splitter. 

Next, we have to include in the consideration the ef- 
fect of the atmosphere. It can be performed using the 
approach of fluctuating- loss channels [llj, [isl - fTsj . Let Og 
and aio be the annihilation operators of the signal and 
the local oscillator, respectively, at the sender. The corre- 
sponding input-output relation for light passing through 
the atmosphere reads as 



bio = Taio + Vl - T^cio, 



(6) 
(7) 



where Cs and Qo are operators of the environment modes 
being in the vacuum state and T is the atmospheric trans- 
mission coefficient. The following properties of relations 
(|ni) and (O are important. First, the transmission co- 
efficient T is a random variable. Second, since the de- 
polarization effect of the atmosphere is negligible, the 
transmission coefficients in both polarization modes are 
perfectly correlated and equal T. Third, in the consid- 
ered case the absence of the depolarization means the 
absence of dephasing, and hence T can be considered as 
a real random variable. Finally, the commutation rules 
require that T e [0, 1]. 

The above treatment can be easily used in Eq. ([T]) with 
the Glauber-Sudarshan P representation [19] for the sig- 
nal density operator. 
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where P{a) is the Glauber-Sudarshan P function of the 
input signal field (at the sender) and |a) is a coherent- 
state vector. Substituting Eqs. ©-dS]) into Eq. ([T]) and 
taking into account that T is a random variable, one gets 
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where 
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(10) 



is the Husimi-Kano Q symbol ^ of the POVM of pho- 
tocount differences, where 

6*1,2 = {r^ + lap ± 2r Re[ae-''^]) , (11) 

V{T) is the probability distribution of the transmission 
coefficient (PDTC) of the atmospheric channel and Iam 



is the modified Bessel function. For simplicity, further 
we refer to the Q symbol of the POVM as the POVM. 

For the purposes of balanced homodyne detection one 
usually uses a local oscillator, which is strong compared 
to the signal. After transmission through the atmo- 
sphere, the intensity of the signal is still small compared 
to the local oscillator, |a|^ <C T^r^. The contribution 
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FIG. 1. (Color online) Scheme of homodyne detection for quantum light passing through the turbulent atmosphere as reported 
in Ref. [ij]. The signal and the local oscillator are sent in two orthogonally polarized modes. After passing through the 
atmospheric channel and collection by a telescope, homodyne detection is realized with a half-wave plate HWP, a polarization 
beam-splitter PBS, and the detectors Di and D2. 



of the latter can be comparable with the noise counts, 
T^r^ ~ Nnc- Following the argumentation of Ref. p^ . 
one can approximate Eq. (jlOp by 
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X exp 
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(12) 



Equations ^ and ((T^ can be directly used for evaluating 
the statistics of photocount differences when both the 
signal and the local oscillator pass through the turbulent 
atmosphere. 



III. POSTPROCESSING NOISE 

The next problem is to connect the photocount differ- 
ences An with the field quadrature, 



(13) 



The corresponding relation is given by (see, e.g., Ref. 0]), 

Afi 



fout 
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where rout is the local-oscillator amplitude scaled by at- 
mosphere losses and detection efficiency. Its real value 
is 



?'out = Ty/rjr. 



(15) 



However, since in the considered experiment we do not 
have any information about the current value of the fluc- 
tuating transmission coefficient T, we can use a certain 
reference value Tjcf , and set 



rout = Trcfy/rjr. 



(16) 



Because the real value of rout given by Eq. (flSl) randomly 
changes in the atmosphere and deviates from its refer- 
ence value pop , the obtained quadrature value suffers 



from a kind of noise, which we refer to as post-processing 
noise. For this reason, the reconstructed density opera- 
tor and any characteristics obtained from the approxi- 
mate quadrature values are, in fact, contaminated by the 
corresponding noise effects. 

Let us consider in more detail the effects of the postpro- 
cessing noise. Based on Eqs. and ([HI), the 
quadrature distribution in the considered case is given by 

+00 

p{x;ip)= / d^aPia) Knoisy[x{f) la] , (17) 



where 
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is the resulting noisy quadrature POVM. Alternatively, 
Eq. (jl7p can be rewritten in the Schrodinger picture as 
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p{x;ip) = / d^aPnoisyia) K[x{(p) -,0] , (19) 



where 

K[x{ip);a] 



: exp 
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is the noiseless quadrature POVM, 
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is the noisy P function of the detected signal, and 
Aa^ gaRoa + d'Hma Laplace operator in phase 
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space. Equation ([21]) can be considered as the quantum- 
state input-output relation, where the noisy density op- 
erator, represented by the P function, is affected by (i) 
fluctuating losses due to the signal transmission through 
the atmosphere, (ii) detection losses and noise counts, 
(iii) postprocessing noise caused by imperfect knowl- 
edge of the transmission coefficient. Any reconstruction 
of the density operator using homodyne-detection data, 
obtained from Eq. (jl4[) together with the approxima- 
tion (|16|) . will yield the noisy quantum state ([21]) . Simi- 
larly, any characteristics obtained from such quadratures 
also correspond to the noisy density operator. 

As follows from Eq. (|2ip . the contribution of noise 
counts in the considered scheme can always be made neg- 
ligible by choosing a sufficiently strong local oscillator. 
In this case, the quantum-state input-output relation re- 
duces to 



noisy 
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This equation can be interpreted as the input-output re- 
lation of a fluctuating- loss channel (cf. Ref. [ll|) with 
the effective transmission coefficient 



Tcs — 



rcf 



(23) 



Additionally, the measurement procedure suffers from a 
kind of effective noise counts whose mean value is 



^rcf 



(24) 



However, the upper bound of Tcff is not restricted any- 
more by the value 1. Similarly, A'off may take negative 
values. As a result, the noisy density operator obtained 
by using Pnoisy(a) in Eq. ([5]) may fail to obey the re- 
quirement of positive semidefiniteness. Such an unusual 
result simply reflects quantum physical inconsistencies of 
the method of data post processing under consideration. 

We consider an illustration of this fact for a single- 
photon-added thermal state (SPATS). This state is ob- 
tained from the single-mode thermal state with the mean 
photon number nth by adding a photon by using para- 
metric down-conversion. The corresponding P function. 
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is regular, which allows its experimental reconstruc- 
tion [2l[. We also assume that fluctuating losses are 
caused by beam wandering only; see the appendix [^ 
This results in log-negative Weibull distribution for the 
PDTC [cf. Eq. ([XT]) ]. 

For a consistent positive-definite density operator, the 
diagonal elements in the Fock-number basis, that is. 



photon-number distribution, should be always nonneg- 
ative. This distribution can be reconstructed from the 
homodyne-detection data; see Refs. and [l^. In the 
strong-turbulence regime, the postprocessing noise may 
result in negative values of the photon-number distribu- 
tion; see Fig. ^a). However, such fake effects can be 
substantially reduced in the case of weak turbulence; see 
Fig.lHb). 
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FIG. 2. (Color online) Photon-number distribution (diagonal 
elements of the effective density operator) for the SPATS with 
nth=l.ll (such as realized in Ref. 21]) disturbed by beam 
wandering for the scheme with unmonitored transmission co- 
efficient. In the corresponding turbulence model, cf. the ap- 
pendixl^ the beam-spot radius, W , on the receiver aperture 
of radius a is 1^=0. 9a, and the standard deviations of beam 
deflection are (a) a=a (strong turbulence) and (b) a— 0.5a 
(weak turbulence). The corresponding reference transmission 
coefficient is T,^t = ^f(T^ f» 0.586 and T,^t = ^/(T^ ^ 0.822, 
respectively. For both cases we use the detection efficiency 
ri—0.5. Negative probabilities demonstrate fake effects caused 
by imperfect postprocessing. 



IV. MONITORED TURBULENCE 

In order to exclude unwanted and fake effects of the 
postprocessing noise, we propose to modify the scheme 
in Fig.dl For this purpose, the signal and the local oscil- 
lator can be split, and a part of the local oscillator can be 
used for monitoring the current value of the transmission 
coefficient; see Fig. O In the most general case the mon- 
itoring will also be affected by different kinds of noise, 
for example by the shot noise of the detector D3. In 
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terms of the previous section, this means that the value 
Tj-ef can now be replaced by the measured transmission 
coefficient Tmcas, which randomly fluctuates and corre- 
lates with fluctuating values of T. The corresponding 
noisy density operator is obtained similarly to Eq. (|2ip. 
However, the PDTC must now be replaced with the joint 
PDTC 'P(r, Tineas) and integrated over both variables T 
and Tineas- The joint PDTC can be given as 

V{T, Trnc^s) - V{T) T'lTncaslT) , (26) 

I 



where 'P(TincaslT') is the probability distribution of the 
measured transmission coefficient under the condition 
that its real value is T. This implies that the input- 
output relation for the considered experimental scheme 
reads as 
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(27) 



where R2 is the reflection coefficient of the beam splitter 
BS2. 

Let us consider the situation when the noise of mon- 
itoring is caused by the shot- noise of the detector D3. 
In this case, the detected number of photocounts is 
related to the measured transmission coefficient Tmcas as 
na = r2ry3|T2pT,^eas + ^3- Here rj^ and are the effi- 
ciency and the mean number of noise counts of detector 
D3, respectively; T2 is the transmission coefficient of the 
beam splitter BS2. The measured transmission coeffi- 
cient Tijioas is thus obtained from the number of photo- 
counts 713 via 



-^mcas (^^3 ) 



1 n,. - TV, 
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(28) 



Since detector D3 records a fraction of the local oscillator, 
rig obeys the Poissonian statistics 



(29) 

This is the probability to get the photocount number 77.3, 
conditioned on the value T of the transmission coefficient. 

The conditional PDTC 'P(TmoaslT) is obtained, by 
transforming random variables 713 — > Tmeas and using 
Eq. (1251), in the form 

-\-oo 

7'(Tnicas|T) = ^ P ("3 1^") ^ [rmcas - Tmcas ("-3 )] • (30) 



From this distribution the conditional expectation value 
E of T,^eas is found to be equal to T^, 

E(t2_|T) = T^. (31) 

The corresponding conditional variance reads as 

Var(Tl_|T) ^ , ^ + ' ^ . (32) 

??3 IT2I H r]^ IT2I r4 



Equations (l^J) and yield the relative error of T^, 
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(33) 



From this expression, it follows that for the measurement 
of T^ with the relative error e one has to use the local- 
oscillator amplitude 



(34) 




+ e27V3. 



Hence a problem of monitoring appears for tiny T val- 
ues, for which the value of the local-oscillator amplitude 
should be really large. However, for this domain the sig- 
nal is of poor quality anyway. We can choose a minimal 
value Tinin of the transmission coefficient T and post- 
select the values T > Tmin. Based on this assumption, 
Eq. (POI) with the relative error e for T^ reduces to 



(^mcas |-^) — ^ (^ -^mcas) : 



(35) 



in which connection the local-oscillator amplitude should 
be chosen according to Eq. (l34| with T=Tijiin. Under 
these conditions the shot noise error of detector D3 be- 
comes small. 

Utilizing the (5-function form psp for the conditional 
PDTC in Eq. one gets 



fnoisy(a)= / dT7'(T) 



(36) 
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The effect of noise counts can be omitted if the local- 
oscillator amplitude obeys the condition 
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(37) 
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FIG. 3. (Color online) Homodyne detection of quantum light passing through the turbulent atmosphere with monitoring the 
transmission coefficient. The signal and the local oscillator are sent in two orthogonally polarized modes. After passing through 
the atmospheric channel and collection by a telescope, the signal and the local oscillator are split by the polarization beam- 
splitter PBS. A part of the local oscillator transmitted through the beam-splitter BS2 is used for monitoring the transmission 
coefficient with the detector D3. Another part, after setting the needed phase by the phase modulator PM and conversion of 
the polarization direction by the half-wave plate IfWP, is combined with the signal on a standard homodyne detector, which 
consists of the 50 : 50 beam splitter BSi and two photodetectors Di and D2. 



In this case, the input-output relation (l36l) reduces to 



(38) 



which appears to be similar to the case of an indepen- 
dently controlled local oscillator, (cf. Ref. ll|). How- 
ever, an important difference is that this relation does not 
contain phase noise of the signal after passing through 
the turbulent channel. It is also worth noting that the S- 
function form of the conditional PDTC ([35l) excludes the 
postprocessing noise of the measured data. For this rea- 
son, no fake quantum effects appear in Eqs. (l36l) and (l38t . 

As in the previous section, we illustrate the 
method with the single-photon-added thermal states [cf. 
Eq. ((25|) ]. Besides, we suppose that this state is displaced 
with the coherent amplitude 7 (cf. Ref. 11]) such that 
its P function is 
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FIG. 4. (Color onhne) P function of the displaced SPATS for 
lma=0, with nth = 1.11 (such as realized in Ref. 21]) and dif- 
ferent values of 7, disturbed by beam wandering and detected 
with monitored transmission coefficient. In the correspond- 
ing turbulence model (cf. the appendix E)) ■ the beam-spot 
radius, W, on the receiver aperture of radius a is W=0.9a, 
and the standard deviation of beam-deflection is a=10a. The 
detection efficiency is rj—O.b. 
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As has been shown in Ref. 'll'], increasing the displace- 
ment amplitude 7 results in diminishing the nonclassi- 
cality in the scenario when the signal and local oscillator 
are radiated from different sources. It turns out that this 
rule does not apply in the considered case. The absence 
of phase fluctuations in input-output relation psp lifts 
strong restrictions for the coherent amplitude of nonclas- 
sical states considered in Ref. [11]. The quantum state 
in this case may preserve its nonclassical properties. In 
Fig. m we show the P function for the scenario with the 
monitored transmission coefficient. It is clearly seen that 
state is still nonclassical even for large values of the co- 
herent displacement amplitude 7. 



V. QUADRATURE SQUEEZING 

Quadrature squeezing is a remarkable property of 
quantum light, which can be observed by homodyne de- 
tection. In this section, we consider how the disturbance 
of the signal and the local oscillator by the atmosphere 
affects the detection of quadrature squeezing. We con- 
sider two orthogonal quadratures [cf. Eq. (fT3|)]. for a 
certain value of the local-oscillator phase (/?, 



Xi = 
^2 = 

A well known relation. 
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(42) 
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i,j=l, 2, connects the covariance matrix with its nor- 
mally ordered form. If a diagonal element (variance) of 
the latter becomes negative for properly chosen ip, the 
state is quadrature squeezed. 

Let us first consider the case with monitored transmis- 
sion coefHcient; see Sec. IIVI Based on Eq. ([55]) . we can 
write the input-output relation for the covariance matrix, 



(43) 



The first term of this relation resembles the standard at- 
tenuation. The second term is caused by the atmospheric 
turbulence. The third term of the equation describes 
the disturbance effect of noise counts on the quadrature 
squeezing. 

The contribution from noise counts disappears for a 
sufficiently strong local oscillator when 



However, due to large contributions of events with small 
T, the value of {T~^) can be really large. For exam- 
ple, if fluctuating losses are caused by beam wandering 
(cf. the appendix 1X1 and Ref. [IBl), this term is infinite. 
In practice, however, the measured transmission coeffi- 
cient is bounded by its minimal value T^i^. This implies 
that the third term of Eq. (|43l) becomes negligible for the 
local-oscillator amplitude satisfying condition ([57)1 . 

It is readily seen from the second term of Eq. (|43l) that 
the disturbance effect of the turbulence on the quadra- 
ture squeezing increases with increasing mean value of 
the quadrature {xi). This means that states with a small 
coherent amplitude have better chances of preserving this 
nonclassical property. For the states with (xi) =0, the 
squeezing of the ith quadrature is disturbed in the same 
way as for standard attenuation. 

Next we consider the case of non-monitored transmis- 
sion coefficient; see Sec. IIIII The corresponding input- 
output relation for the covariance matrix is obtained from 
Eq. (I^T]) and reads as 



r, (T^s) (: Ax.Ax, :) + {x,) (x,) (AT^^) 

/ (reff)-l iVnc 

I 2 + r;r2r2^ 



(44) 



Si J, 



where Tcs is given by Eq. ((23)) . In contrast to the case 
of monitored transmission coefficient, the third term in 
this equation does not disappear even in the case of weak 
noise counts. 

The post-processing noise may result in such an effec- 
tive covariance matrix that the corresponding density op- 
erator is not positive semidefinite anymore. This means 
that the Heisenberg uncertainty relation. 



is not satisfied; see Fig. [5] When the real transmission 
coefficient appears to be much less compared with the ref- 
erence transmission coefficient Tj-ef , the negative effective 
noise [cf. Eq. (PH)] leads to such a nonphysical squeezing. 
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FIG. 5. (Color online) Product of noisy variances, 
(Aa;?) . (Axy) . , vs the reference transmission coefB- 

\ ' noisy \ ' noisy 

cient, Tref, for the scheme with non-monitored turbulence. 
The solid (dashed) line corresponds to the vacuum state (8- 
dB squeezed vacuum state) at the transmitter. We suppose 
that fluctuating losses are caused by beam wandering with 
the standard deviation of beam deflection cr=40a and with 
the beam-spot radius W— 0.95a, which leads to 35-dB mean 
losses. The detection efficiency is 77=0.5. The shaded area 
corresponds to the violation of the Heisenberg uncertainty re- 
lation. 

The covariance matrix completely characterizes quan- 
tum states in the Gaussian approximation. They can 
be used, for example, for continuous-variable protocols 
of quantum-key distribution with coherent and squeezed 
states [2^. Evidently, in order to get a consistent effec- 
tive density operator in this approximation, one should 
use an appropriate reference transmission coefficient Tj-ef . 
This means that the effect of the post-processing noise 
may preserve consistency of the noisy density operator 
even in the case of using the scheme with non-monitored 
transmission coefficient. 
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VI. SUMMARY AND CONCLUSIONS 

Homodyne detection of quantum light passing through 
the turbulent atmosphere with the local oscillator sent 
in the orthogonally polarized mode is a promising tech- 
nique for long-distance quantum communication based 
on continuous variables. In this connection, a problem 
appears of how to connect a measured photocount differ- 
ence with the field quadrature. Indeed, when the receiv- 
ing local-oscillator amplitude is a fluctuating variable, 
this question is not trivial anymore. We consider two 
possible solutions for this problem. One possibility could 
be based on some reference value for the local-oscillator 
amplitude transmission coefficient (e.g., its mean value). 

(45) 

Alternatively, here we proposed a method based on the 
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monitoring of the transmitted local oscillator with the 
aim of having control over the fluctuating transmission 
coefficient of the turbulent atmosphere. In both cases, 
the quantum state of the received light can be charac- 
terized by a noisy density operator, which includes also 
information about shortcomings of the measurement and 
postprocessing procedures of the used methods. 

When the local-oscillator amplitude (and thus the 
transmission coefficient) is monitored, the main limita- 
tions are caused by stray-light and dark-count noise. 
These effects can be, in principle, eliminated by a suf- 
ficiently strong local oscillator and the postselection of 
events with an appropriately chosen threshold value of 
the transmission coefficient. In the simpler procedure, 
using a fixed reference value of the local-oscillator am- 
plitude, the shortcomings caused by the resulting post- 
processing noise are much more dramatic. The resulting 
disadvantages cannot be eliminated anymore. In such 
a scenario, the noisy density operator may even violate 
the fundamental requirement of positive semidefiniteness. 
Thus fake quantum effects may occur due to the used 
post-processing procedure. However, even based on this 
technique one may obtain a consistent noisy density op- 
erator in the Gaussian approximation, provided that the 
reference transmission coefficient is properly chosen. We 
believe that these methods may be of some interest in the 
context of continuous- variable quantum key distribution. 
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ating losses are caused by beam wandering. If the beam is 
randomly deflected around the aperture center according 
to a two-dimensional Gaussian distribution with the vari- 
ance cr^, the PDTC is given by the log- negative Weibull 
distribution, 



(Al) 



for T e [0, To] and else. Here the parameters Tq, A, and 
R are expressed in terms of the beam-spot radius at the 
aperture W and the aperture radius a. 



A =5 



lo = V 1 - exp 



exp[-4^]li(4^. 



In 



(r 



2Ti 



exp[-4^]Io(4f, 



R 



In 



1 



exp[-4^]Io(4^) 



(A2) 



(A3) 



(A4) 



In the case when the turbulence is weak and the beam 
is focused on the aperture, the beam-deflection variance 
can be approximately evaluated as 
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Appendix A: PDTC for beam wandering 

Here we remind readers of some results of Ref . (16i | , in 
particular the explicit form of the PDTC when the fluctu- 



1.919 Clz^i2Wor^^^, 



(A5) 



where is the index-of-refraction structure constant. 
Wo is the beam-spot radius at the radiation source, 
and z is the distance between source and receiver aper- 
ture 'F, '2?]. Integration with the PDTC V{T) must be 
performed in the Lebesgue sense with respect to the mea- 



sure d 



i?(21n^); 
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